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Positive vacuum energy together with extra dimensions of space imply that our four- 
dimensional Universe is unstable, generically to decompactification of the extra dimen- 
sions. Either quantum tunneling or thermal fluctuations carry one past a barrier into the 
decompactifying regime. We give an overview of this process, and examine the subsequent 
expansion into the higher-dimensional geometry. This is governed by certain fixed-point 
solutions of the evolution equations, which are studied for both positive and negative 
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we also outline a possible mechanism for compactification to a four-dimensional de Sitter 
cosmology. 
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1. Introduction 



We now have good reason to believe that we live in an accelerating Universe; this 
point has particularly been brought home with the recent WMAP results 0, combined 
with earlier cosmological observations. It is also widely believed that our fundamental 
description of nature should involve extra small dimensions of space. These two statements 
alone lead one to a very general argument that the Universe as we know it is unstable to 
a catastrophic transition^]. 

The generic instability is for the extra dimensions of space to begin to grow, and 
our world to evolve into a higher-dimensional one. However, depending on the details 
of the potential for the shape and size moduli of the extra dimensions, there may also 
be basins of attraction with negative potential, which lead to equally catastrophic big 
crunches. The instability towards expansion of the extra dimensions was first argued 
for in the context of string theory by Dine and Seiberg||. Their arguments were based 
on supersymmetry. However, || points out that the underlying mechanism is a simple 
dynamical one, driven by the dynamics of long-distance gravity, and is independent of the 
existence of supersymmetry. 

Examples of potentials exhibiting this instability are now being widely studied in 
string theory. Flat directions in moduli space have been a long-standing problem in string 
theory. Recent developments in compactifications with fluxes and branes[f4[-|l4| have pro- 
vided examples of dynamics that fixes these moduli [ l^JTG ] . In particular, |16|] showed that 



by adding an anti-D3 brane to the solutions of ||15|| , one can lift the vacuum energy and 
find locally stable minima with positive cosmological constant. Subsequent work on other 
approaches to vacua with positive cosmological constant has included [ ]TT| - |iTi[ . 

These developments have led to much discussion, as well as criticism[|0|,^l|, of the 
resulting picture of a "landscape" of stringy vacua. This picture has led to a forceful res- 
urrection [[22] of the idea that constants of nature - particularly the cosmological constant - 
are determined anthropically; the large number of possible fluxes and resulting vacua |23| , p4 | 



(for a review see |[25|| ) together with the observation that in a sufficiently large distribution 



one expects to find a small enough cosmological constant have given strong fuel to this 
possibility. 



One of the criticisms from |20| is potentially relevant to our discussion, and so deserves 
comment. Banks and Dine argue that an effective potential description of the landscape 
is not justified and cannot be trusted. Underlying this argument is the realization that 
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there are regions in the combined field space of the moduli and metrics where the dynamics 
becomes strongly coupled. While it is certainly not inconceivable that this could render the 
entire picture inconsistent, we take a more sanguine perspective. It may well be that there 
are dragons off in the mountains of the landscape. However, the valley we find ourselves 
in seems perfectly tame, and we expect that dynamics nearby is likewise tame. Of course 
we'd very much like to understand the strongly coupled dynamics to understand how we 
arrived to our present vacuum, but for now we'll take the perspective that the dynamics of 
the full quantum wavefunction has somehow deposited us here, and our problem is to see 
what happens next. We provisionally accept that the effective potential is a useful tool in 
this investigation. 

In the next section we review the derivation of the effective potential for the radial 
dilaton modulus. In particular, this analysis shows that modular landscapes have a generic 
feature, much like the "front range" of Colorado - the mountains taper off to a semi-infinite 
plain. As we roll into this region, the extra dimensions of space expand. Section two also 
discusses examples of potentials that can be obtained from fluxes and branes, for example 
in string theory. Section three turns to the problem of analyzing the asymptotic dynam- 
ics of solutions that have escaped a metastable minimum and are running to infinitely- 
expanded extra dimensions. Asymptotically these solutions become fixed-point solutions, 
whose form we derive from the equations of motion. Section four discusses mechanisms for 
escape from the metastable minimum - thermal excitation over or tunneling through the 
barrier. Finally, in section five, we assemble these results with a general discussion of the 
decompactification transition out of (our?) metastable de Sitter space and into the decom- 
pactifying regime, with the resulting higher-dimensional evolution. We also briefly discuss 
a scenario in which transitions can take place back and forth between a higher-dimensional 
dS space and the four-dimensional one; in this case, the dominant configuration should be 
determined by the relative entropies, providing a possible mechanism for compact ificat ion. 
An appendix provides a more detailed phase-plane analysis, for both positive and negative 
spatial curvatures, of the solutions arising in our discussion of decompactification, 

Spontaneous decompactification has also been recently discussed for simple compact- 
ifications with fluxes in [p7| |, which appeared while our paper was being written up. 
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2. Radial dilaton dynamics: generalities 

2.1. Dimensional reduction 

We begin by discussing dynamics of the radial dilaton, following^]. Specifically, sup- 
pose that we begin with (ci+4)-dimensional action 

S = J d d+4 xV^G M d+2 TZ + C{ip) + t{i),K) , (2.1) 

where X and G are the coordinates and metric of the full (<i+4)-dimensional spacetime, 
Mp is the (<i+4)-dimensional Planck mass, TZ is the Ricci scalar, C(ijj) is the lagrangian 
representing the leading contribution of generic matter sources in a derivative expansion, 
possibly including localized sources such as D-branes, and C(ip, TZ) summarizes possible 
corrections to the leading lagrangian that involve higher powers of the curvature and/or 
higher derivatives of matter fields. This action may be the effective action for string theory, 
or for some other fundamental theory of gravity. For our present purposes we assume that 
all the moduli except the overall volume modulus of the internal space are fixed, e.g., as 
in ||15|| ; we will consider the coupled dynamics of this modulus and the four-dimensional 



metric. (It has been argued in |28j that in the case of multiple moduli with exponential 



potentials, a single term dominates, leading to an obvious generalization of our analysis to 
the multi-moduli case.) Thus we assume that ( |2.1|) has solutions of the form@ 

ds 2 = ds 2 4 + R 2 (x)g mn (y)dy rn dy n . (2.2) 

The action governing solutions with R{x) = e D ^ varying slowly on the compactifi- 



cation scale follows from dimensional reduction of (|2.1| ). The Einstein-Hilbert term gives 
a 4d effective action 

S EH = M d+2 V d J d^x^l [e dD(x) ^ 4 + d(d - 1) ( VD) 2 e dD{x) + e (d " 2)D ^ d ] , (2.3) 

where Vd and IZd are the volume and curvature of the <i-dimensional compact metric g mn . 
However, in this frame the 4d effective Planck mass varies with R. To make contact with 
usual treatments of 4d dynamics, one should choose a new set of units by performing the 
Weyl rescaling 

9^u -> e- dD g 4 ^ ; (2.4) 



1 Here we suppress a possible warp factor, which should not change our general picture. 
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the 4d Planck mass then becomes 



Ml = Mp 2 V d . (2.5) 

The additional terms in (|2.1|) then contribute an additional potential for D. The net result 
is an action of the form 

S = Ml J d A x^~ A \n A - ^d(d + 2)(VD) 2 -V(D)^j . (2.6) 

A central observation of || is that any physics that stabilizes the radial dilaton from 
runaway to infinite volume, D — > oo, must do so only locally - the Weyl rescaling ( |2.4| ) 
introduces an inverse power of the volume-squared into the potential. This means that 
to stabilize the radial dilaton, the higher-dimensional dynamics would have to produce an 
energy density growing at least as fast as the internal volume for large volume. This does 
not appear realizable in any realistic physical theory. 

V 
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Fig. 2: A sketch of a potential with a metastable de Sitter region, and a 
minimum with negative cosmological constant. If the Universe fluctuates into 
the AdS basin of attraction, it evolves to a big crunch singularity. 

Therefore if there is a positive potential minimum representing the present de Sit- 
ter phase of our Universe, this minimum must be metastable, and the potential should 
generically appear as in fig. 1 or fig. 2. Fig. 2 produces an instability to a big crunch 
spacetime |]29| , j30! . Our focus will be the more generic case of runaway to infinitely ex- 
tended extra dimensions, as follows from fig. 1. 

Of course, in full generality there will be a nontrivial multi-dimensional moduli space of 
compact manifolds, and a potential function on this space. This has recently been explored 
in the case of string theory for the region of string configuration space corresponding to 



flux compactifications, and the resulting configuration space dubbed [p2| 'the landscape.' 
The argument given in 0] implies that the landscape has a general feature similar to the 
geology of the front range of the Rocky Mountains in North America: the mountains and 
valleys of the landscape roll off into a flat plain, extending to infinity. 

2.2. The dilaton potential - examples 

In general the dilaton potential receives contributions from non-trivial field config- 
urations in the extra dimensions. For illustration, it's particularly fruitful to consider 
examples provided from string theory, particularly branes and fluxes. There may well also 
be other effects, both perturbative and non-perturbative in the couplings. 
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Some of the examples from string theory were enumerated in [0. A g-form flux field 
with action 

r f 2 

S p (x- J d d+4 XV^G-± , (2.7) 

gives a potential 

V F oc R~ d ~ 2q (2.8) 
and a space-filling p-brane with action 

S p = / dV p+1 , (2.9) 

9s Jm 4 xC p - 3 

gives a potential 

V p oc RP- 3 ~ 2d . (2.10) 



From the action ( |2.3|) , curvature on the internal manifold can give a potential 

V R oc R- d ~ 2 . (2.11) 



Nonperturbative effects give other contributions, e.g., [|T6| included a potential of the form 

V = Be~ 2aRA /R s , (2.12) 
with different constants a and s arising from euclidean D3 brane instantons or gluino 



condensation. Finally a' corrections were argued in |3T[| to produce potential contributions 
of the form 

Vk oc . (2.13) 



3. Asymptotic dilaton dynamics 

3.1. Fixed-point dynamics 

Quantitative details of the de Sitter decay will depend on the contribution of these 
and other possible terms to the potential. However, we'll be particularly interested in the 
asymptotic structure of the expanding region of higher-dimensional space. The dynamics 
inside this region is governed by the asymptotics of the potential at large dilaton, and for 
this in general we need only to focus on the leading term. Since this will typically be of 
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the form V oc R n , let us analyze more closely the dilaton dynamics with a single such term 
in the potential. With canonically normalized dilaton, 



<f)= ^d(d + 2)D , (3.1) 



we find an action of the form 

S = M't I d*Xy/=g~ 

with 



^4-^W) 2 -^) 



(3.2) 



V&) = \<T a + ■ (3.3) 



In these conventions, 



a = ' l dT2\ 1 + 2 i ] q - flux 



d ( 1+ d + l-P ) p _ brane (3.4) 



d + 2 

d + 2 
d 



internal curvature 



Friedman-Robertson- Walker solutions for actions of this form, with exponential po- 
tentials, have been investigated previously in the literature, see e.g., [|32| - p9| . Recently, 
there has been a great deal of discussion of analogous solutions in string- and M-theory, 
see e.g., PD| - |5B| , P7| . As we'll see in the next section, such FRW solutions, with spatial 
curvature k = 0, ±1, are useful in studying the decompactification transition. These have 
the form 

ds\ = -dt 2 + a 2 (t)dx 2 , <f> = <f>(t) (3.5) 

where d\ 2 is the appropriate spatial metric with curvature k = 0, ±1. The equations of 
motion take the form 

a 2 6 12 6 



2 „+7^ + -> (3-6) 



Here we've included the possibility of a stress tensor with a barotropic equation of state, 



Pj = (7 - l )Pi ) (3-9) 

for which the density evolves as 

p\ = — • (3.10) 

a 

In the case k = ±1, the dominant effective contribution to p at large a is that summarizing 
the spatial curvature, which gives 

P2/3 = ~% ■ (3.H) 

Massive matter (7 = 1) or radiation (7 = 4/3) will make subdominant contributions as 
a — > 00 in these cases, but if k = their contribution can be relevant to this asymptotic 
behavior. 

The asymptotic behavior of solutions of these equations is typically governed by fixed 
point "tracker" solutions. The fixed points relevant to the cases k = 0, — 1 have been found 
via phase space analyses in p5| , |37|1 . A simple extension of this work covers the case k = +1 
as well, as described in the appendix. Let's understand how these solutions arise. The 
starting point is the observation that for a scale factor that increases as a power of time 
(which will be justified shortly), 

a = ao(ct) p , (3.12) 



and for an exponential potential (|3.3| ), equation (|3.8j ) drives to infinity as a power: 

= {ctf (3.13) 

where the constant c is fixed to be 



c=-J-^. (3.14) 
2 y 3/3 - 1 K 1 

Note that ( |3.8| ) can also be integrated to give the ^dependent terms in the Freedman 
equation: 



2 V = J_ 
12 ' 6 ~ a 2 t 

Thus, using ( |3.10|) , the Freedman equation (|3.6| ) becomes 

d 2 c' , 



h + ^ = ^- (3-15) 



a 2 t^P a 2 t 2 



+ ^ (3-16) 



for some constant d . First consider the cases k = 0, — 1, where d is positive. Which term 
gives the dominant fixed-point asymptotic behavior as t — > oo depends on the relative 
magnitudes of a and 7. A phase plane analysis of the dynamics appears in the appendix, 
but the basic features can be understood straightforwardly directly from the equations. 
For 37/9 > 2, the first term becomes subdominant as t — > 00. Then, using the power law 
form of the solution, (|3.12|) , we find f3 = 1/a 2 . This means the condition for this fixed 



point to dominate is a 2 < 3^/2. In contrast, for a 2 > 3^/2 both terms are relevant. 
Correspondingly, the power law form ( |3.12| ) now implies (3 = 2/37. The complete phase 



plane analysis (see [|35| , |37H and the appendix) shows that these are indeed the fixed points 
to which the generic solution is attracted — in the appendix, fixed point (a) for a 2 < 3^/2 
and (b) for a 2 > 3^/2 (see fig. 4). Any other fixed points are either unstable nodes or 
saddle points. 

The case of positive spatial curvature k = +1 allows negative d and is somewhat 
more subtle. For a 2 < 3^/2 = 1, depending on initial conditions, the solution may reach 
a turning point and reflect to a collapsing universe, a < 0, rather than reaching the 
asymptotic a — > 00 region. If this happens, other terms in the Freedman equation ( |3.6| ), 
e.g., due to matter/radiation fields, become important in the dynamics, and a general 
analysis is not possible. However, if the solution asymptotically expands, the analysis is 
similar to the k = — 1 case. For a 2 > 1, solutions generically reach a turning point and 
recollapse; this behavior is not well-described by the power law form ( |3.12| ) . These features 
are well illustrated by the corresponding phase plot in the appendix — see fig. 5. 

We can summarize the resulting asymptotically expanding solutions by 

ds 2 = -dt 2 + a 2 G [c(t - t )] 2(3 dx 2 , (3.17) 

and 

e* = [c(t-t )] 2/a , (3.18) 



with c given in (|3.14|). The constant ao may be fixed by equation (|3.6|) for k = ±1 and the 



expansion exponent (3 is given by 



/3 = max( A, _L) . (3.19) 



3.2. Specific cases 



Consider first the case of nonzero spatial curvature, k = ±1. At large radius, the cur- 
vature dominates over the matter contributions to the Friedmann equations, and curvature 
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corresponds to 7 = 2/3. Thus for a > 1, which is generically the case for string-theory 
induced potentials, for k = — 1 we have (3 = 1, giving a non-accelerating solution with 
a linearly growing scale factoril; for k = +1 we have collapse. For example, consider a 
flux-induced potential. We can easily see from ( |3T4] ) that 



<8»=(l + 2i)( 1+ 2|^] (3.20) 



which gives a 2 > 1. Similarly for internal curvature, we find 

2 



^curvature — 1 + ^ (3-21) 



so again a 2 > 1. The case of a brane-induced potential gives, from (|3[ 



For p < <i + 2, z.e., a brane of codimension one or higher, we therefore also find a 2 > 1. 
An exceptional case is that of a completely space-filling brane, p = d + 3, or equivalently 
higher-dimensional cosmological constant, which gives 



a =VlT2 <1 - < 3 ' 23) 

The corresponding solution then has (3 = 1/a 2 = 1 + 2/c? and is accelerating. 

Other string theory potentials typically fall off even more rapidly at infinity, removing 
them further from this accelerating case. 

Next consider spatially flat universes, k = 0. In these cases, the asymptotic density 
of matter or radiation, and in particular the dominant value of 7, is relevant for deter- 
mining which attractor governs the dynamics. Note, however, that one still only achieves 
acceleration for a < 1. 

We will discuss the .D-dimensional interpretation of these solutions in the following 
section. 



However, solutions approaching this fixed point may be eternally accelerating 1 40 > 48 , 50 1 
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4. The decompactification transition 

4-1. Escaping inflation 

We next turn to the study of decay of a metastable de Sitter minimum in the potential 
for the radial dilaton. In particular, we assume that the relevant dynamics for our discus- 
sion is that of the radial modulus, and that other possible moduli of the internal manifold 
are fixed, although as we've pointed out this analysis should extend to the multi-moduli 
case where that is dominated by an exponential potential due to a single modulus, as in 
||28|| . Thus we consider the lagrangian (|3.2| ), where V(<p) has a form sketched in fig. 1. For 
simplicity, we work in units where M4 = 1 in the following. 

Any given region of the de Sitter universe corresponding to the metastable minimum 
will ultimately decay to a decompactifying universe. The time scale depends on the pa- 
rameters of the potential.i If the minimum corresponds to our presently inflating phase, 
we should have Vq = V(<po) ~ 10 -120 <C 1. For a general potential, without any further 
fine tuning, we expect that the height of maximum of the barrier at <p± is V\ ~ 1 and that 
the width of the barrier is also A0 ~ 1 . There are then two basic mechanisms for inducing 
spontaneous decompactification. The first is for thermal excitations of de Sitter space to 
take one over the top of the barrier. The second is quantum tunneling through the barrier, 
as described by Coleman and de Luccia [|9| . 



The thermal activation rate is given by the action of the Hawking-Moss instanton||54T 



as has been argued from the stochastic approach to inflation [pD] - pT| . This gives the proba- 
bility of a thermal fluctuation to take one to the top of the barrier. If the change in entropy 
for this fluctuation is AS, the probability is given by the usual formula P = exp{AS}. 
This probability is 

^thermal = e**>>- S <*> , (4.1) 

where where S(cpo) = — 247r 2 /V"o is the action of the de Sitter instanton, ie, the four- 
sphere solution, with <p = <po, the metastable minimum. Similarly, S(<pi) is the solution 
with (p = (pi, the maximum of the barrier. The de Sitter recurrence time|6^] is given by 
T r = ex.p{S (<po)} . The thermal decay lifetime is thus shorter than this by an exponential 
factor: 

^thermal = 6 ^ 1 T r . (4-2) 



For a more detailed review of parameters and dynamics see |16|. 
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The tunneling probability is given by [29.16f 



P tU nnel = eX P | [i + (4V3r2)]2 j (4.3) 

where the tension of the bubble wall is 

T= / d<P^2V(<P) . (4.4) 



For the parameters that we generically expect, we see that Vq <T 2 . The probability 
therefore yields a decay time 

e -64^/T 2 Tr . (4.5) 



Comparing ( [4.2| ) and ( |4.5|) shows that which process is dominant clearly depends on the 
relative magnitudes of the tension T and the barrier height V\. 

Once the field excites past the barrier, it classically evolves towards decompactifica- 
tion. In the case of thermal activation, one expects a picture where an entire horizon 
volume thermally excites over the barrier. The boundary conditions then start the solu- 
tion near the maximum of the potential, and the resulting classical solution can roll into 
the decompactification region. However, fluctuations are expected to be important in this 
case, and so one cannot clearly identify the resulting dynamics as a k = 0, ±1 universe.^ 
However, these may serve as rough indicators for the dynamics, bearing in mind that 
different regions may have different spatial curvatures. One also expects to have some 
excited matter fields in the resulting solution; this is important to avoid certain collapsing 
solutions in the case of vanishing spatial curvature, k = 0, but more generally it seems 
plausible that collapse could occur in some regions, and expansion in others. 

In the case of tunneling, one can be a little more specific about the resulting boundary 
conditions, by matching the CdL instanton onto a lorentzian geometry at the turning point. 
The 0(4) symmetry of the instanton continues to an SO(3, 1) symmetry of the subsequent 
classical solution, implying [f29j|63| that one has evolution on spatial sections with k = — 1. 



These features are illustrated in fig. 3. The boundary conditions on the "initial slice" of 
the cosmology are 

(j) = , a = , d = l and <fi = 4>t (4-6) 
where 4> t is the value at the turning point. 



We thank to A. Linde for a discussion on this point. 
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Fig. 3: A representation of the the Coleman de-Luccia instanton, describing 
tunneling from a metastable de Sitter region into a bubble of space in which 
the extra dimensions of space are expanding. The lower half of the diagram is 



the euclidean solution of [29], for the potential V((f>). This matches onto the 
lorentzian solution pictured in the upper half of the diagram. The straight 
lines correspond to the a = surface of the resulting expanding k = — 1 cos- 
mology; the surfaces above these are surfaces of constant t in this cosmology. 
The boundary of this decompactifying region is a bubble wall (circular below, 
hyperbolic above) which asymptotically expands at the speed of light into the 
metastable de Sitter region. 

4-2. Expansion to higher dimensions 

In either case, even if we know the potential explicitly it is typically hard to find an 
explicit solution. However, the fixed-point behavior of these systems tells us that this is not 
necessary to understand the asymptotic behavior. Specifically, we expect the asymptotic 
behavior of the decompactifying solution to be governed by the leading term in the potential 
as a — > oo. This is expected to generically be one of the exponentials we've considered, 
and so decompactification is asymptotically described by one of the fixed-point solutions 
of the preceding section, described by (in the k = 0, — 1 cases) the solution (|3.17 )- (|3".19|) . 



Since, as shown by eqn. (|3.18|) , the compact manifold expands, the relevant description 



becomes the higher-dimensional description rather than that of four dimensions. Recall 
that these differ by the rescaling (|2.4| ), so in terms of the fundamental units of the higher- 
dimensional theory, the asymptotic metric takes the form 

ds 2 4+d = e- dD ® {-dt 2 + a 2 (t)d X 2 ) + e 2D ^ds 2 . (4.7) 
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In discussing the higher-dimensional form of these solutions, it proves easier to work 
with the coefficient A, related to a by 



\ d + 2 
d a 

thus for fluxes, branes and internal curvature, we find from equation (|3.4| ) 

2q 



A = 1 + 



d 



q — flux 



= l + d+3 - p p-bra„e 



d 



(4.9) 



= 1+ d 



internal curvature . 



Then given the asymptotic solution in equations ( 3.17 ) and (|3.18| ), the higher-dimensional 
metric ( |4.7[ ) becomes 

ds 2 d+4 = (c(t - t ))- 2/X [-dt 2 + a 2 (c(t - t )) 2p d X 2 ] + (c(t - t )) 4/Xd ds 2 . (4.10) 
For A > 1, this metric can be simplified by defining a new time coordinate 



1 A 



r = - 



c A- 1 



(c(t-V)) 



l-l/A 



(4.11) 



with which we find 



ds = —dr + a (cr) 



Here the constant c is given by 



2 jl „2 fA s2(/3\-l)/(X-l) ,2 



dx 2 + (&rf d ^ ds 2 . 



A- 1 



(4.12) 



(4.13) 



recall that c was given in eq. (|3.14| ), and f3 in eq. (|3.19|) . Both the three large dimensions 
and the compact dimensions expand. However, the compact dimensions expand more 
slowly unless 

(/?A-l)d<2. (4.14) 

Thus for k = — 1 they do so in the case of flux-induced potentials and of potentials arising 
from branes of codimension greater than two. 

The special case of a wrapped p = d + 3 space-filling brane corresponds to A = 1 and 
(3 = 1 + 2/d. Now we define 

(4.15) 



r = - ln(c(t - t )) 
c 



14 



and the metric becomes 



ds 2 = -dr 2 + e AcT ' d (a 2 d X 2 + ds 2 d ) . (4.16) 

Thus, locally, the metric looks like a patch of inflating deSitter space — for a discussion 
of such asymmetric foliations of deSitter space, see . 

As the size of the compact dimensions become comparable to that of the other three, 
in general the four-dimensional effective description will break down and one must de- 
scribe the solutions as fully higher-dimensional solutions. Since the field configurations, 
e.g., branes and fluxes, on the compact manifold will generically be no n- uniform, this 
means that the resulting higher-dimensional solution is generically non-uniform. More- 
over, fluctuations of modes that get light in this limit can become relevant. Nonetheless, 
in many cases we expect the four-dimensional solutions that we have described to lift to 
give a higher-dimensional solution, and thus a reasonably accurate picture of the higher- 
dimensional dynamics. 



5. Overview: the fall 

Let us assemble the pieces of the picture that we have discussed. Beginning with 
our Universe in its currently inflating phase, ultimately a fluctuation will carry it out of 
our metastable minimum. This may either be a thermal fluctuation over the barrier, or 
a tunneling event through the barrier. In the thermal case, an entire horizon-sized region 
fluctuates over the barrier. In the tunneling bubble forms, and then expands. 

The relevant timescales are extremely long, as they contain a factor of the recurrence 
time, T r ~ exp{10 120 }. The decay time, expressed as a fraction of the recurrence time, is 
given in eqns. (|4.2|) , Q4.5Q , and depends on the relative magnitude of the bubble tension 
and the barrier height. 

What happens next depends on which basin of attraction the fluctuation takes us into. 
There may be accessible basins of attraction with negative effective cosmological constant; 
in that case the Universe will undergo a big crunch. Alternatively, we know that there 
is generically an infinite basin of attraction corresponding to decompactification. This 
is the one that we have considered in this paper. We also expect this basic picture to 



extend to the case of multi-moduli runaway; ref. |28| has argued that the case of multiple 
runaway scalars with exponential potentials can be reduced to the single-scalar case. There 
also may be basins of attraction where some of the dimensions of space decompactify and 
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others remain stabilized. These would be described by straightforward generalization of 
our analysis. 

Once we find ourselves in such a basin of attraction, the Universe decompactifies, 
with asymptotic dynamics typically given by (|4.12| ) or (|4.16|) . In the case of a thermal 
fluctuation, we expect an entire horizon-sized region to evolve into the region where higher- 
dimensional dynamics is relevant, but the actual configuration may be highly non-uniform, 
with collapse in some regions and expansion in others. It is in this context that we imagine 
that our results for k=+l might be applicable. That is, they might describe the dynamics 
of a small patch with positive spatial curvature in a larger inhomogeneous solution. The 
surprising result from the phase plane analysis - - see appendix - - is that for typical 
potentials, e.g., induced by fluxes or branes, the solutions generically evolve to a big crunch. 
We expect that such regions collapse into black holes and are thus casually disconnected 
from the expanding regions.! 

In the case of tunneling, a bubble of higher-dimensional space forms, and its walls 
expand into the four-dimensional region at the speed of light. The solution inside the 
bubble is given by the asymptotic dynamics that we have described. As in old inflation, 



the majority of space continues to be inflating four-dimensional space |6^,^5[, but any given 
point will ultimately transition into the higher-dimensional realm. 

Now we return to a special case of of space-filling brane with p = d+3, i.e., a (positive) 
cosmological constant in (d+4) dimensions. In this case, the landscape as illustrated in fig. 
1 is deceptive. Despite the appearance of the usual infinite plain, there are no solutions 
rolling towards flat decompactified space in (d+4) dimensions. Rather as we have explicitly 
shown in ( [116|) , the (apparently) metastable dS space transitions to a spacetime which 
asymptotes to (<i+4)-dimensional de Sitter space. Recall from the discussion around ( |3.23| ) 
that this was the special case for which the universe appeared to be accelerating from a 
4d point of view. 

Hence this introduces the following interesting possibility. Here we may consider the 
transition our four-dimensional dS space to a higher-dimensional dS space. However, in 
this case, we also expect the reverse process to be possible. Moreover, we expect that, if 
the system can indeed be thought of as a thermal ensemble and the de Sitter entropies as 
accurate representations of the number of states, the solution that dominates the ensemble 



A similar description has been advocated for the crunch singularities appearing in transitions 



towards AdS minima [22]. 
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should be that with the higher entropy. This could be either the higher-dimensional or 
lower-dimensional de Sitter space depending on the relative magnitudes of the entropies 

M (d+4)(d+2)/2 

Sd + 4 P Ad+2)/2 ( 5 - 1 ) 

and 

S4 ~ a7 X4 — ' (5 ' 2) 

where Va is the volume of the compact manifold at the "metastable" minimum. In fact, 
it is not difficult to construct examples where the compactified solution has the larger 
entropy In particular, we know that a minimum corresponding to the present universe 



has an entropy S4 ~ 10 120 , and so if the hypothetical higher-dimensional theory had a 
cosmological constant that is order unity in Planck units, our four-dimensional metastable 
minimum would be expected to strongly dominate such an ensemble. Hence this scenario 
could provide a mechanism for compactification from the higher-dimensional theory. 

Certainly, this is not the entire story. Typically within this framework, compactified 



anti-de Sitter solutions will also arise, e.g., by adjusting fluxes[]66[|. It is not clear what 
their role should be in such a scenario [p0| , |21|j30f l ■ However, this aspect of the dynamics 
seems a fruitful ground for future research. 
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Appendix A. Phase plane analysis for k = ±1 

Following the phase plane analysis of [[371 , we describe the general solutions of the 



equations (]3.6|) - (|3.8|) for the cases where the spatial metric is either positively or negatively 
curved. Hence in equations ( |3T6| ) and ( |3.7| ), we have 7 = 2/3 and substitute for p 1 as in 
equation (|3.10| ). The phase plane dynamics explicitly reveals the presence of fixed point 
'tracker' solutions, discussed in the main text, as the generic end-points of various solutions. 
We begin by defining 



A\JSa/a voa/a a 

With this choice of variables, the Freedman constraint ( |3.6| ) becomes 

x 2 + y 2 -kz 2 = l (A.2) 



which can be used to eliminate z. The second-order equations (|377|) and ( p^q ) can be 
written in the form: 

x' = -2x(l - x 2 ) + y 2 (VSa - x) , (A.3) 

y' = y (1 - y 2 ) - xy (73a - 2x) , (A.4) 

where prime denotes d/d log a and a is the (positive) constant characterizing the exponen- 
tial potential, as in equation (|3.3|). Note that these equations are independent of k. The 
feature which then distinguishes these two cases is the range of the variables as determined 
by the constraint (|A.2 ), i.e., x 2 + y 2 < 1 for k = — 1 and x 2 + y 2 > 1 for k = +1. We 



are also generally interested in solutions describing an expanding universe and so we focus 
our attention on y > 0. Since eqns. ( |A.3|) , ( |A4j) are symmetric under y — > —y, the y < 
dynamics is a simple extrapolation. 

These equations exhibit the following fixed points: 

a. x = x a = a/y/Z, y = y a = a/1 - « 2 A z 2 = 0, 

b. x = Xb = l/\/3a, y = yb = V^/VSa, z 2 = k (a~ 2 — l), 

c. x = 0, y = 0, z 2 = —k, 

d. x = 1, y = 0, z 2 = 0, 

e. x = —1, y = 0, z 2 = 0. 

Which of these fixed points is physically realized (and their stability properties) depends 
on the values of parameters k and a. In particular, the fixed point (a) only appears for 
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Fig. 4: Fixed points and some typical trajectories for k = — 1 in the regimes: 
(i) a < 1, (ii) 1 < a < \/3, and (iii) a > y/3. (The specific plots above were 
made for a = 1/V3, 5/2y3, and 2.) 
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a < a/3 in order that y a is real. While (b) has Xb and y^ real for all values of a, it is only 
relevant for a < 1 with k = +1 and a > 1 with /c = — 1 so that the fixed point lies in 
the appropriate domain of the phase plane, i.e., z 2 > 0. Similarly, (c) is only relevant for 
k = —1. Finally the fixed points (d) and (e) always appear relevant independent of the 
parameters. Figs. 4,5 illustrate the fixed points and the flows in the phase plane for the 
various distinct parameter ranges. We now comment briefly on these results. 

Consider the case k = — 1, as described by Fig. 4. First for a < 1, (a) is the only 
stable fixed point and then describes the end-point behavior of generic solutions. As 
described in the text, amongst the string theory potentials, this situation only seems to 
be realized with a completely space-filling brane (with p = d + 3), in which case this fixed 
point corresponds to an asymptotically de Sitter solution in the full (d+4) -dimensional 
spacetime. For 1 < a < a/3, (a) becomes a saddle-point and then disappears (i.e., moves 
off into the complex plane) for a > v^3- Hence in the regime a > 1 which seems to be the 
generic case for string theory, (b) is the stable fixed point.i 

Fig. 5 illustrates the flows for k = +1. In this case, (a) is again a stable fixed point 
for a < 1 and the same comments given above apply here. In fact, (b) only appears in this 
parameter range as well, but it is an unstable saddle-point. Hence for a > 1, the generic 
flows run off to infinity in directions bounded by y 2 /x 2 < 2. This behavior is also generic 
for a broad range of initial conditions when a < 1. Physically these flows correspond to 
expanding cosmologies where the kinetic energy in the scalar field dominates the potential, 
and the spatial curvature leads to decelerated expansion. As a result a = 1/z reaches zero 
at some finite value of a (and t). Hence beyond this point, they begin collapsing. In 
the phase plane, this transition is realized by mapping the asymptotic flow with (x, y) to 
(—a;, —y) and following the flow backwards with respect to the arrows indicated in the 
figures. In the present context, the solutions of interest in the present context are those 
which run off to infinity with x > 0, i.e., a, > 0. In this case, for any value of a, they 
flow towards a big crunch at the fixed point (e). However, as commented in the main text, 
we expect other terms, e.g., matter and perturbations, to become important in the final 
stage of this crunch. 

Note that during the collapse stage above, we have d, x < and hence <j) > 0. That 
is, the internal space continues to expand for these solutions. From a higher-dimensional 

6 Qualitatively, these results apply for the case of any barotropic fluid with a positive energy 
density [j37f — the latter would be the dominant energy contribution for the case k = 0. 
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Fig. 5: Fixed points and some typical trajectories for k = +1 in the regimes: 
(i) a < 1, (ii) 1 < a < \/3, and (iii) a > \f?>. (The specific plots above were 
made for a = 1/2, 5/2^3, and 2.) 
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perspective, however, we would still regard this as a collapsing phase, as given the full 
metric Q4.7| ), one can easily show that the proper volume element is contracting both in 
the noncompact four-dimensions and in the full (<i+4)-dimensional solution. 

We have focussed on the end-point of the decompactification trajectories, but it is 
interesting to consider the full trajectories and follow the flows backwards to an initial 
fix-point. Generically, for the entire range of parameters, the flows originate at one of the 
repulsive fixed points (d) or (e). Hence these solutions emerge from a singular big bang. 
As with the crunch above, the analysis provided here is expected to be incomplete as many 
other terms will be important near this initial singularity. The generic appearance of these 
singularities was, of course, observed by [pO| , pT|j30[| and used as part of their criticism of 
string theory landscape picture. However, for tunneling from a metastable dS vacuum, 
these singularities don't appear as part of the evolution of the system. Another possible 
viewpoint on this may be to consider an "extended" landscape, in which the 4d radius (or 
even more general metric degrees of freedom) is included; in this space the singular regions 
may be well separated from the solutions of interest, e.g., by being in different basins of 
attraction, though this question deserves further exploration. 

We might also comment that, in any range of parameters considered above, there are 
many trajectories which begin with negative x and cross over to positive x. These would 



fall into the class of M-theory solutions, which were recently discussed |4jj - |49[] because they 
exhibit an accelerating phase. As first noted in ||42j| . near x = 0, the scalar potential 
momentarily dominates the energy density and so cosmic acceleration is seen by four- 
dimensional observers. However, this is only a momentary phase since, as noted in the 
main text, generically the M-theory potentials are two steep to generate an (continuously) 
accelerating solution. 
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